Inherent stochasticity of superconductive-resistive switching in nanowires 
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Hysteresis in the current-voltage characteristic 
in a superconducting nanowire reflects an under- 
lying bistability. As the current is ramped up 
repeatedly, the state switches from a supercon- 
ductive to a resistive one, doing so at random cur- 
rent values below the equilibrium critical current. 
Can a single phase-slip event somewhere along 
the wire — during which the order-parameter fluc- 
tuates to zero — induce such switching, via the lo- 
cal heating it causes? We address this and related 
issues by constructing a stochastic model for the 
time-evolution of the temperature in a nanowire 
whose ends are maintained at a fixed tempera- 
ture. The model indicates that although, in gen- 
eral, several phase-slip events are necessary to in- 
duce switching, there is indeed a temperature- 
and current-range for which a single event is suf- 
ficient. It also indicates that the statistical dis- 
tribution of switching currents initially broadens, 
as the temperature is reduced. Only at lower 
temperatures does this distribution show the nar- 
rowing with cooling naively expected for resis- 
tive fluctuations consisting of phase slips that are 
thermally activated. 

The essential qualitative characteristics of quasi-one- 
dimensional superconducting nanowires are controlled 
by fluctuations of the superconducting order parameter, 
these fluctuations being predominantly thermal or quan- 
tal, depending on the temperature regime^. Bulk su- 
perconductors undergo a sharp transition from an elec- 
trically resistanceless (i.e. superconducting) to a resis- 
tive (i.e. normal) state, e.g., with increasing tempera- 
ture. In contrast, as explained by Little^ and Langer 
and Ambegaokar 3 , in quasi-one-dimensional supercon- 
ductors the resistanceless (and truly long-range-ordered) 
state is destabilized by a certain class of accessible order- 
parameter fluctuations that connect topologically dis- 
tinct sectors of current-carrying states. In so doing, these 
fluctuations, which are known as phase-slip events, can 
dissipate supercurrent, and because of them such sys- 
tems undergo a broad evolution between the (nominally) 
superconductive and normal states, e.g., with increasing 
temperature. 

Recent advances in the fabrication of ultra-narrow su- 
perconducting wires — using carbon nanotube— or DNA- 
templating^ — have spurred renewed interest in one- 
dimensional superconductivity and opened up new av- 
enues for investigating the impact of order-parameter 
fluctuations. One setting in which order-parameter fluc- 
tuations in superconducting nanowires have been widely 
investigated, both theoretically and experimentally, is 
that of transport properties in the vicinity of the normal- 
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FIG. 1: Model a. Schematic of an experimental configura- 
tion described by our model: a superconducting nanowire is 
suspended between two thermal baths, b. Sketch showing the 
attenuation of the order parameter in the core of a phase-slip, 
c. Schematic of the simplified model. All phase slips are taken 
to occur in the central (i.e. shaded) segment of length I, which 
is assumed to be at a uniform temperature T; heat is carried 
away through the end segments, which are assumed to have 
no heat capacity. The temperature at the ends of the wire is 
fixed to be Tt,. d. Sketch of a typical temperature profile. 



to-superconducting quantum phase transition^ - — . In 
this setting, the primary mechanism underlying destruc- 
tion of (nominal) superconducting order is depairing as- 
sociated with magnetic fields or magnetic impurities. As 
is well known, applied currents also cause depairing and, 
if larger than a certain value (known as the thermody- 
namic critical or depairing current), would render the 
superconducting state locally unstable (regardless of the 
role of phase-slip fluctuations}^ - —. 

However, phase-slip fluctuations, which are responsible 
for the broad resistive transition in quasi-one-dimensional 
superconductors, also allow for premature switchingi 3 - - — 
to the resistive state, i.e. a nonequilibrium transition 
from the (nominally) superconducting, low-resistivity 
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FIG. 2: Hysteresis in the I-V characteristic obtained by find- 
ing the steady-state solutions of Eq. JTJ in the continuous 
Joule heating limit. 

state to the (nominally) normal, high-resistivity one. If 
damping of the order-parameter dynamics were low, a 
single phase-slip event would induce such switching, in 
analogy with what happens in undcrdamped Josephson 
junctions. By contrast, nanowires are generally over- 
damped, and so, whilst causing resistance, phase slip- 
page does not, by itself, induce switching. As discussed 
in Ref. (T^Ji this resistance causes Joule heating which, 
if not overcome sufficiently rapidly by conductive cool- 
ing, effectively reduces the depairing current, ultimately 
to below the applied current, thus causing switching to 
the highly resistive state. Naturally, this switching is 
not deterministic, owing to the underlying stochasticity 
of the phase-slip events that are responsible for the re- 
sistance. Rather, for a given subcritical applied current 
there is a statistical distribution of times at which switch- 
ing occurs, characterized by a mean switching time (i.e. a 
superconducting state "lifetime"). 

Our focus here is on stochastic aspects of the 
superconducting-to-resistive switching dynamics, an area 
that has not received much attention, to date. Inter 
alia, by obtaining the current-dependent mean switch- 
ing time and convolving it with the sweep rate of the 
applied current that describes the experimental proto- 
col, we shall determine the statistical distribution of 
currents at which switching occurs. Besides its fun- 
damental significance, the characterization of switching 
dynamics in nanowires seems likely to have technolog- 
ical implications, such as for the integration of super- 
conducting wires into electronic circuitry as controllable 
(current-limiting) switching elements, the implementa- 
tion of nanowire-based devices ^ 16 ' 17 , and the exploration 
of the use of nanowires in quantum computers. 

Having in mind the configuration in recent and ongoing 
experiments on superconducting nanowires, we consider 
a free-standing wire of length L and cross-sectional area 
A, the ends of which are held at a fixed temperature 
lb, as shown in Fig. [1] The fact that the wire is free- 
standing (i.e. lacks any substrate) is conducive to a clear 



interpretation of the measurements. On the other hand, 
the absence of an overall thermal bath means that any 
heat generated locally in the wire by a source term Q can 
be taken away only through the ends; the corresponding 
heat conduction equation for the temperature 0(x,i) at 
position x along the wire at time t reads 

C v (Q)d t e(x,t)=d x [K s (e)d x Q(x,t)]+Q(x,t), (1) 

and is characterized by the specific heat C v (0) and ther- 
mal conductivity K s (<d) of the wire, together with the 
boundary condition 0(±L/2,t) = 7b at its ends. Note 
that although our analysis rests on the premise that there 
are no additional heat-removing channels, it can readily 
be extended to account for such possibilities. 

Before addressing dynamical issues, let us dwell briefly 
on the steady-state solutions of the heat conduction equa- 
tion, obtained by setting d t <d(x,t) = and assuming 
that the wire is subjected to temporally continuous Joule 
heating at a rate given by ALQ(x) — I 2 R(Q(x), I). Here, 
the function R(Q',I) is to be understood as the resis- 
tance of an entire wire held at a uniform temperature 
Q(x) = 0'. The system- wide I-V characteristic at a 
given boundary value It, can be traced by obtaining the 
temperature profile 0(x) for every value of current I in 
both up and down (parametric) sweeps of /. On deter- 
mining i?(0, I) via the current-biased version of LAMH 
theor y 3 ' 18 i 19 for phase-slips, the I-V curves are indeed 
found to become progressively more hysteretic in I as Tb 
is lowered (see Fig. [2]). This steady-state problem was 
previously studied by Tinkham et alii*, who used for 
i?(0, 1) the experimental linear-response resistance mea- 
sured at Tb = O, leading to qualitative agreement with 
the hysteresis observed in MoGe nanowires^. 

Our aim here is to study the inherent stochasticity in 
the switching process, and therefore it is necessary for us 
to explicitly take into consideration the fact that the re- 
sistive fluctuations of the superconducting nanowire con- 
sist of discrete phase-slip events (labelled by i) that take 
place at random moments of time U and are centered at 
random spatial locations Xi. The work done on the wire 
by a phase slip may be obtained from the time integral 
of IV (t), in which the Josephson relation d<fi/dt = 2eV/h 
may be used to relate the voltage pulse to the rate of 
change of the phase difference^, via fundamental con- 
stants H and e. Hence, a single phase slip (or anti-phase 
slip), which corresponds to a decrease (or increase) of <fi 
by 27T, will heat (or cool) the wire by a "quantum" of en- 
ergy hi /2e. Thus we arrive at the central thrust of our 
paper: the dynamics of switching from the superconduct- 
ing to the resistive state in the nanowires is controlled by 
a heat conduction equation that is stochastic by virtue 
of its source term: 

Q(x,t) = ^^^(TiFix-x^t-U), (2) 

i 

where F(x — x{) is a normalized (to unity) form factor 
representing the relative spatial distribution of heat pro- 
duced by the i th phase-slip event, and <7j = ±1 for phase 
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(anti-phase) slips. The probability per unit time T± for 
an anti-phase (phase) slip to take place, depends on the 
local temperature 0(x,i) and the current /. The ran- 
domness in Xi and ti generates the stochasticity in the 
switching from the superconductive to the resistive state. 

To capture the essential physics whilst making the 
problem amenable to analysis, we shall consider the sim- 
pler model, represented in Figs. [Tfc and [Hi. Given that 
the edge effects favor phase-slip locations away from the 
wire ends, the source term is restricted to the region 
near the center of the wire. The system is thus mod- 
eled by assuming that (i) the heating takes place within 
a central segment of length I to which a uniform tem- 
perature T is assigned, and (ii) the heat is conducted 
away through the end segments, within which we ignore 
the heat capacity^. To simplify the problem further, we 
make use of the fact that T + <C T_ and ignore the pro- 
cess of cooling by anti-phase slips. To account indirectly 
for their presence, we use a reduced rate T = T_ — T + 
instead of T_ for phase-slip events. This ensures that the 
discrete expression for Q will correctly reduce to the con- 
tinuous Joule-heating expression, in view of the LAMH 
formula R(T,I) = hT/2eI. 

By using the model defined above, the description re- 
duces to a stochastic ordinary differential equation for 
the time-evolution of the temperature of the central seg- 
ment: 



dT 
dt 



-a(T, T b )(T - T b ) + r)(T, I) £ 6(t - U), (3) 



where the second term on the RHS corresponds to heat- 
ing by phase slips, and the first term to cooling as a re- 
sult of conduction of heat from the central segment to the 
external bath via the two end-segments, each of length 
(L — I) /2. The temperature-dependent cooling rate a is 
given by 



a(T,T b ) 



l(L - l)C v (T) T — 



l — dT'K 3 {T>). (4) 



If Ti and Tf are temperatures before and after a phase 
slip then, using 



r Tt hi 
A I J t C v (T')cff'=|i, 



(5) 



we can express the temperature 'impulse' due to a phase 
slip,i.e Tf — Ti = ^(T, I) = rj(T{, I), as function of either 
Ti or Tf, depending on the context. 

Let us now elucidate the physical and mathematical 
structure of Eq. ([3]). To begin with, we shall consider the 
continuous-heating limit, rj{T 1 1)T(T, I), for the source 
term, and express Eq. ([3]) as dT/dt = —dU/dT. In 
Fig.[3j we illustrate the form of the 'potential' U(T, T b , /) 
for fixed T b : there is a range of currents / for which U 
has two local minima, corresponding to the supercon- 
ducting (at low-T) and the resistive (at high-T) states, 
separated by a local maximum. The resulting bistability 
is central to the underlying physics. On the one hand, 
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FIG. 3: Effective potential U(T,T b ,I) (dashed line) and the 
mean first-passage time r as functions of the temperature 
T of the central segment for various bias currents I and for 
Tb = 1.2 K. The marks on the temperature axis indicate 
the temperatures that the central segment would have after 
1,2, 10 phase slips in the absence of cooling (as given by 
Eq. © for To = T b ). 



it explains the origin of the hysteretic behavior; on the 
other hand, it provides a basis for phrasing the ques- 
tion of stochastic switching dynamics in superconducting 
nanowires in terms of an existing general framework for 
stochastic bistable systems. In what follows, we focus on 
the stochastic variable T(t); to ease the notation we do 
not display the dependences on / and T b unless essential. 
To continue the analysis of the stochastic equation, imag- 
ine turning off the cooling term. If we now start with an 
initial temperature Tq then 



T , T + tKTj), To + 77(T ) + TyfTo + V (T )), 



(6) 



defines the discrete sequence of values that T jumps to, 
as marked on the horizontal axes in Fig. [3] for To = T b . 
The probability per unit time, T(T), to make a jump 
changes at each step, and so does the size rj(T) of the 
jump, owing to their explicit dependence on temperature. 
On the other hand, if we turn off the heating term then 
we have a deterministic problem in which T would decay 
at a rate a(T), from its initial value To > T b to the 
bath temperature T b , which is the lowest value T can 
have. It is the competition between the discrete heating 
and the continuous cooling that makes for a rather rich 
stochastic problem. We hope that our solution will also 
furnish insight into other physical problems that possess 
a similar mathematical structure. 

The master equation for P(T,t), the probability for 
the temperature of the central segment of the nanowire 
to be T at a time t (given that it had some initial value 
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To at time to), reads 

d t P(T, t) = d T [(T - T b ) a(T) P(T, t)] - T(T) P(T, t) 
+ T(T - rj(T)) P(T~ rj(T), t)(l-cfr rj(T)) , 



(7) 



where the first (i.e. the transport) term corresponds to 
the effect of cooling, and the last two terms correspond 
to the effects of heating. Note that the term (1 — drrjiT)) 
appears because of the dependence of the jump size on T, 
as given by rj(T) . The fundamental quantity of interest 
is the mean switching time T s (Tb,I), i.e. the mean time 
required for the wire to switch from being superconduc- 
tive to resistive, assuming that the wire has temperature 
T = Tb when the current I is turned on at time t = 0. 
The master equation, Eq. (|7|) , provides the starting point 
for generalizing the standard procedure for computing r s 
via the evaluation of the mean first-passage time^I. 

The mean first-passage time r(T — > T*), to go past a 
point T = T for the first time having started from some 
T < T , can be shown to satisfy the equation 

-(T h -T)a(T)d T T(T)+T(T) [t(T) - t(T + r,(T))] = 1, 

, (8) 

together with the conditions r(T) = for T > T and 
dr(T)/dT = at T = Tb, which are appropriate for our 
problem. Some illustrative plots for r(T — > T ), obtained 
by numerically solving Eq. ([5]) are shown in Fig. [3J with 
the choice of T being somewhat larger than the location 
of the local maximum of U. From these plots we see 
that the mean first-passage time has a plateau at low 
values of T and then rapidly decreases in the vicinity 
of the potential barrier. At very high currents, as can 
be seen from the last panel of Fig. [3j the local stability 
of the superconducting state disappears, and so does the 
plateau in the mean first-passage time. In these plots, the 
tick marks on the T axes correspond to the temperatures 
given by the sequence ((6]) for To = 2b. 

As long as the high-T minimum is lower than the 
low-T one, and T is chosen to be appreciably past 
the intervening potential maximum (in order to elim- 
inate the possibility of reversion to the superconduct- 
ing state), we can make the identification r s (Tb,I) = 
r(Tb — > T , Tb,I). The number of tick marks (see se- 
quence ([6])) between Tb and T is nothing but the of 
number N(Tb,I) of phase-slip events required to raise 
the temperature of the central segment from Tb to T 
in the absence of cooling. Accordingly, N(Tt,,I) also 
provides an estimate of the number of phase-slip events 
needed to overcome the potential barrier if the timespan 
of these events is insufficient to allow significant cooling 
to occur. 'Thermal runaway' — heating by rare sequences 
of closely-spaced phase slips that overcome the potential 
barrier — constitutes the mechanism of superconductive- 
to-resistive switching within our model. As the N(Tb, I) 
becomes large, the total number of phase-slip events tak- 
ing place before switching can happen, and correspond- 
ingly the value of r s (Tb,I), may indeed be quite large. 
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FIG. 4: Switching statistics a. Map showing N(Th,I) (see 
text) and the contour lines (solid lines) for the inverse of 
mean switching time, r s _1 = 1, 10 3 , 10 6 s _1 ; the contour lines 
(dashed lines) for the phase-slip rate F are also shown. The 
thermodynamic (depairing) critical current (dashed-dotted 
line) is plotted for reference, b. Switching-current distri- 
butions Psw obtained at various values of Tb and for r = 
58/uA/s. c. The logarithms of r s _1 (colored lines) and of V 
(thinner black lines) as a function of /, obtained for the same 
set of T> values as in panel (a). The colors of r s _1 plots cor- 
respond to different values of N(T^,I) [as indicated in the 
legend of panel (a)]. 



Our key findings are summarized in Fig. [3J There is a 
region of / and Tb for which the occurrence of just one 
phase slip is sufficient to cause the nanowire to switch 
from the superconductive to the resistive stated; in this 
case TjT 1 = r. A switching measurement in this range 
can thus provide a way of detecting and probing a sin- 
gle phase-slip fluctuation. As, outside this range, several 
phase-slip events are required for switching, r s _1 deviates 
from r (see panel A graphical representation of the 
contour lines for a few values of rr 1 and T, chosen in an 
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experimentally accessible range, is provided in panel [4^. 
Whilst the spacing between the T contour lines decreases 
monotonically on lowering Tb, the spacing between r s _1 
lines can be seen to behave non-monotonically. 

The mean switching time r s in bistable current-biased 
systems can be either directly measured or extracted 
from the switching-current statistics^ generated via the 
repeated tracing of the I-V characteristic by ramping 
the current up and down at some sweep rate r. For this 
reason, in Fig. HJd we have illustrated the behavior of 
this distribution of switching currents in superconduct- 
ing nanowires based on the theory presented here. Upon 
raising Tb, one would naively expect the distribution to 
become broader for a model involving thermally activated 
phase slips. Such an broadening in the distribution- width 
is indeed obtained up to a crossover temperature scale 
T b =r (r) (i.e. the temperature below which, loosely speak- 
ing, switching is induced by single phase slips). How- 
ever, on continuing to raise Tb, but now through tem- 
peratures above T b 3r (r), the distribution- width shows a 



seemingly anomalous decrease. This is a manifestation of 
the now-decreasing spacing between the r s contour lines. 
This striking behavior above T" may be understood by 
the following reasoning: the larger the number of phase- 
slips in the sequence inducing the superconductive-to- 
resistivc thermal runaway, the smaller the stochasticity 
in the switching process and, hence, the sharper the dis- 
tribution of switching currents. This non-monotonicity in 
the temperature dependence of the width of switching- 
current distribution, along with the existence of a regime 
in which a single phase-slip event can be probed, are the 
two key predictions of our theory 
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